A-type Antiferromagnetic and C-type Orbital-Ordered State in LaMnOa 
Using Cooperative Jahn- Teller Phonons 
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The effect of Jahn- Teller phonons on the magnetic and orbital structure of LaMnC>3 is investigated 
using a combination of relaxation and Monte Carlo techniques on three-dimensional clusters of Mn06 
octahedra. In the physically relevant region of parameter space for LaMnOa, and after including 
small corrections due to tilting effects, the A-type antiferromagnetic and C-type orbital structures 
were stabilized, in agreement with experiments. 
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The theoretical understanding of manganese oxides is 
among the most challenging .-current areas of research 
in condensed matter physics.™ Experimental studies of 
manganites have revealed a rich phase diagram origi- 
nating from the competition between charge, spin, and 
orbital degrees of freedoms. A simple starting frame- 
work for Mn-oxide investigations is contained in the 
double-exchange (DE) ideas, where ferromagnetism in- 
duced by hole doping arises from the optimization of the 
hole kinetic energyfl In addition, recent results using the 
one-orbital model revealed a more complicated ground 
state, with phase separation tendencies strongly compet- 
ing with ferromagnetism,B leading to a potentiaLexpla- 
nation of the Colossal Magneto- Resistance effect J3 

However, to understand the fine details of the phase 
diagram of manganites the one-orbital model is not suf- 
ficient since the highly nontrivial A-type spin antiferro 
(AF) and C-type orbital structures observed experimen- 
tally in the undoped material LaMnOsu cannot be prop- 
erly addressed in such a simple context. Certainly two- 
orbital models are needed to consider the nontrivial state 
of undoped manganites. In this framework the two-band 
model without phonons has been studied before, and the 
importance of the strong Coulomb repulsion has been 
remarked for the appearance of the A-AF stateBt] How- 
ever, Coulombic approaches have presented conflicting 
results regarding the orbital order that coexists with the 
A-type spin state, with several approaches predicting G- 
type orbital order, which is not observed in practice. 
In addition, many experiments suggest that Jahn- Teller 
(JT) phonons are important in manganites and, thus, 
at present it is unclear whether the dominant interac- 
tion between electrons in Mn-oxides should be consid- 
ered Coulombic or phonon mediated. For this reason it 
is important to analyze if a purely JT-phononic calcula- 
tion is able to reproduce the experimental properties of 
undoped manganites, goal that provides the motivation 
for the present paper. The main result observed in this 
effort is that A-AF order, in combination with a C-type 
orbital arrangement, is indeed induced by JT- phonons 
in realistic parameter regions for LaMn03, namely large 
Hund-coupling between e g -electrons and t2 g -spins, small 



AF interaction between t2 g -spins, and strong electron- 
lattice coupling. This shows that JT-based calculations 
can lead to correct qualitative predictions for mangan- 
ites, actually improving on purely Coulombic approaches 
in the orbital sector. 

To carry out the calculations note that experiments 
have revealed an orbital structure tightly related to the 
JT-distortion of the Mn06 octahedron. If each JT- 
distortion would occur independently, optimal orbitals 
can be determined by minimizing the kinetic and in- 
teraction energy of the e g -electrons, as in models with 
only Coulomb interactions. However, oxygens are shared 
between adjacent MnOg octahedra, indicating that the 
JT-distortions occurs cooperatively. Particularly in the 
undoped situation, all MnOe octahedra exhibit JT- 
distortionSj— indicating that such a cooperative effect is 
important. E2I Thus, in order to understand the magnetic 
and orbital structures of LaMnOs, the electron and lat- 
tice systems must be optimized simultaneously. How- 
ever, not much effort has been devoted, to the microscopic 
treatment of the cooperative effect,Eil although the JT- 
effect ip-ithf Mn-oxides has been addressed by several 
groups .EStd To remedy this situation, here a computa- 
tional investigation of cooperative JT-phonons in man- 
ganites is carried out, focusing on the n = 1 density, 
where n is the electron number per site. 

The motion of e g -electrons tightly coupled to the lo- 
calized t2 g -spins and the local distortions of the MnC<6 
octahedra is described by 



+ J' £ s i • Sj + (1/2) YSP<& + Qli + Qli), (i) 

(iJ) i 

where c\ aa [ciba) is the annihilation operator for an e g - 
electron with spin a in the d x 2_ y 2 (d 3z 2_ r 2) orbital at 
site i. The vector connecting nearest-neighbor (NN) sites 
is a, fIL, is the hopping amplitude between 7- and 7'- 
orbitals connecting NN-sites along the a-direction via 
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the oxygen 2p-orbital, Jh is the Hund coupling, Si the 
localized classical t2 g -spin normalized to |Si| = 1, and 
tj = (ci, 02,(73) are the Pauli matrices. The dimension- 
less electron-phonon coupling constant is A, Qu denotes 
the dimcnsionlcss distortion for the breathing mode of 
the MnOg octahedron, Qi\ and are, respectively, JT 
distortions for the (x 2 — y 2 )- and (3z 2 — r 2 )-type modes, 
and <7o is the unit matrix. J' is the AF-coupling between 
NN t2 g -spins, and (3 is a parameter to be defined below. 

To account for the cooperative nature of the JT- 
phonons, the normal coordinates for distortions of the 
MnOg octahedron are writterO as Qu = (l/v / 3)(^xi + 
Lyi + L zi ), Q 2i = (l/V2)(L xi - Lyi), and Q 3i = 
(l/v6)(2L z i — Lxi — L y i), where L a i denotes the distance 
between neighboring oxygens along the a-direction, given 
by L a i = L a + (uf — wf_ a ). Here, L & is the length between 
Mn-ions along the a-axis and uf denotes the deviation of 
oxygen from the equilibrium position along the Mn-Mn 
bond in the a-direction.Ej In general, L a can be different 
for each direction, depending on the bulk properties of 
the lattice. Since the present work focuses on the mi- 
croscopic mechanism for A-AF formation in LaMnOs, 
the undistorted lattice with L x = L y = L z is treated 
first, and then corrections will be added. In the cubic 
undistorted lattice, the hopping amplitudes are given by 
^aa = — V^ab = — V^ba = 3i* b = t for the x-direction, 
^aa = V^^ab = V^t^ a = ^bb = * f° r tne y-direction, 
and tl h = At/3 with t\ & = t* h = = for the z- 
direction. The energy unit is t. The parameter f3 is 
defined as (3 = (wbr/wjiO 2 , where Wbr and wjt are the vi- 
bration energies for manganite breathing- and JT- modes, 
respectively, assuming that the reduced masses for those 
modes are equal. Using experimental results and band- 
calculation data for Wbr and wjt,cj it can be shown that 
(3 w 2. However, the results presented here are basically 
unchanged as long as (3 is larger than unity. 

To study Hamiltonian Eq. (1), two numerical tech- 
niques were here applied. One is the relaxation tech- 
nique, in which the optimal positions of the oxygens are 
determined by minimizing the total energy. In this cal- 
culation, only the stretching mode for the octahedron, 
namely uf + uf_ a = 0, is taken into account. Moreover, 
the relaxation has been performed for fixed structures 
of the t2 g -spins such as ferro (F), A-type AF (A-AF), 
C-type AF (C-AF), and G-type AF (G-AF), shown in 
Fig. 1(a). The advantage of this method is that the op- 
timal orbital structure can be rapidly obtained on small 
clusters. However, the assumptions involved in the relax- 
ation procedure should be checked with an independent 
method. Such a check is performed with_thc unbiased 
MC simulations used before by our groupo The domi- 
nant magnetic and orbital structures are deduced from 
correlation functions. In the MC method, the clusters 
currently reachable are 2x2x2,4x4x2, and 4x4x4. 
In spite of this size limitation, arising from the large num- 
ber of degrees of freedom in the problem, the available 
clusters are sufficient for our mostly qualitative purposes. 



In addition, the remarkable agreement between MC and 
relaxation methods lead us to believe that our results are 
representative of the bulk limit. 




Ferro A-AF C-AF G-AF 




FIG. 1. (a) Magnetic structures in 2 x 2 x 2 clusters, (b) 
Total energy as a function of J' on a 2 x 2 x 2 lattice with 
Jh = 8 and A = 1.5. The solid lines and circles indicate 
the relaxation and MC results, respectively. MC simulations 
have been performed at temperature 1/200. (c) Spin correla- 
tion function S(q) obtained by MC simulations as a function 
of J', at Jh = 8 and A = 1.5. Solid and open symbols denote 
the results in 4 x 4 x 2 and 4x4x4 clusters, respectively. 
Circles, squares, diamonds, and triangles indicates S(q) for 
q = (0, 0, 0), (w, 0, 0), (it, it, 0), and (w, 1T, 7r), respectively, (d) 
Optimized orbital structure for each magnetic structure. 

In Fig. 1(b), the mean-energy is presented as a func- 
tion of J' for J H = 8 and A = 1.5, on a 2 x 2 x 2 clus- 
ter with open boundary conditions. The solid lines and 
symbols indicate the results obtained with the relaxation 
technique and MC simulations, respectively. The agree- 
ment is excellent, showing that the relaxation method 
is accurate. The small deviations between the results of 
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the two techniques are caused by temperature effects. As 
intuitively expected, with increasing J' the optimal mag- 
netic structure changes from ferro- to antiferromagnetic, 
and this occurs in the order F^A-AF^C-AF-^G-AF. To 
check size effects, the t2 g -spin correlation function S(q) 
was calculated also in 4 x 4 x 2 and 4x4x4 clusters, 
where 5(q) = (1/N) £y e-^-j^SrSj), N is the num- 
ber of sites, and (• • •) indicates the thermal average value. 
As shown in Fig. 1(c), with increasing J' the dominant 
correlation changes in the order of q = (0, 0, 0), (ir, 0, 0), 
(tt, 7t, 0), and (jr, it, 7t). The values of J' at which the spin 
structures changes agree well with those in Fig. 1(b). 

The shape of the occupied orbital arrangement with 
the lowest energy for each magnetic structure is in 
Fig. 1(d). For the F-case, the G-type orbital structure 
is naively expected, but actually a more complicated or- 
bital structure is stabilized, indicating the importance of 
the cooperative treatment for JT-phonons. For the A-AF 
state,E3 only the C-type structure is depicted, but the G- 
type structure, obtained by a 7r/2-rotation of the upper 
x-y plane of the C-type state, was found to have exactly 
the same energy. Small corrections will remove this de- 
generacy in favor of the C-type as described below. For 
C- and G-AF, the obtained orbital structures are G- and 
C-types, respectively. Although the same change of the 
magnetic structure due to J' was already reported in the 
electronic model with purely Coulomb interactions,!^ the 
orbital structures in those previous calculations were G-, 
G-, A-, and A-type for the F-, A-AF, C-AF, and G-AF 
spin states, respectively. Note that for the A-AF state, 
of relevance for the undoped manganites, the G-type or- 
der was obtained,EI although in another treatment for the 
Coulomb interaction, the C- and G-type structures were 
found to be degeneratejj as in our calculation. 




FIG. 2. Magnetic phase diagram on the (J', A) plane for (a) 
Jh = 4 and (b) 8 (relaxation method). Below the thin solid 
lines in the F and A-AF regions, the JT-distortion disappears, 
suggesting that the system becomes metallic. 

In Figs. 2 (a) and (b), the phase diagrams on the 
(J', A)-plane are shown for Jh = 4 and 8, respectively. 
The curves are drawn by the relaxation method. As ex- 
pected, the F-region becomes wider with increasing Jjj. 
When A is increased at fixed Jh, the magnetic structure 
changes from F— > A-AF^C-AF^G-AF. This tendency is 
qualitatively understood if the two-site problem is consid- 



ered in the limit Jh 3> 1 and Ejt 3> 1, where Ejt is the 
static JT-energy given by E JT = A 2 /2. The energy-gain 
due to the second-order hopping process of e g -electrons is 
roughly SEaf ~ 1/Jh and SEp ~ 1/-Ejt for AF- and F- 
spin pairs, respectively. Increasing -Ejt, SEp decreases, 
indicating the relative stabilization of the AF-phase. In 
our phase diagram, the A-AF phase appears for A > 1.1 
and J' < 0.15. This region does not depend much on 
Jh, as long as Jh S> 1. Although A seems to be large, 
it is realistic from an experimental viKamoint: Ejt is 
0.25eV from photoemission experimentscJ and t is esti- 
mated as 0.2 ~ 0.5eV,li3 leading to InSnA < 1.6. As for 
J', it is estimated as 0.02 < J' < 0.1.113 Thus, the loca- 
tion in parameter-space of the A-AF state found here is 
reasonable when compared with experimental results for 
LaMn0 3 . 




FIG. 3. (a) Total energy as a function of on the 2x2x2 
lattice with L x = L y > L z for Jh = 8, A = 1.6 and J' = 0.05 
(relaxation method). The solid and dashed curves denote the 
C- and G-orbital states, respectively, (b) Orbital structure 



in the A-AF phase for t% 



and L x — L y > L z . These 



shapes are very close to purely 3x 



and 3y — r types. 



Let us now focus on the orbital structure in the A-AF 
phase. In the cubic lattice studied thus far, the C- and 
G-type orbital structures are degenerate, and it is unclear 
whether the orbital pattern in the x-y plane corresponds 



to the alternation of 3* 



and 3y 2 — r 2 orbitals ob- 



served in experiments.^ To remedy the situation, some 
empirical facts observed in manganites become impor- 
tant: (i) The MnOg octahedra are slightly tilted from 
each other, leading to modifications in the hopping ma- 
trix. Among these modifications, the generation of a 
non-zero value for is important, (ii) The lattice is 
not cubic, but the relation L x w L y > L z holds. From 
experimental resultsp these numbers are estimated as 
L x = L y = 4.12A and L z = 3. 92 A, indicating that the 
distortion with (^-symmetry occurs spontaneously. Note 
that the hopping amplitude and J' along the z-axis are 
different from those in the x-y plane due to this distor- 
tion. Motivated by these observations, the energies for C- 
and G-type orbital structures were recalculated including 
this time a nonzero value for t| a in the magnetic A-AF 
state (see Fig. 3(a)). In the real material, it can be shown 
based on symmetry considerations that the tilting of the 
MnOg octahedra will always lead to a positive value for 
i| a . Then, the results of Fig. 3(a) suggest that the C-type 
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orbital structure should be stabilized in the real materi- 
als, and the explicit shape of the occupied orbitals is 
shown in Fig. 3(b). The experimentally relevant C-type 
structure with the approximate alternation of 3x 2 — r 2 
and 3y 2 — r 2 orbitals is indeed successfully obtained by 
this procedure. Although the octahedron tilting actually 
leads to a change of all hopping amplitudes, effect not 
including in this work, the present analysis is sufficient 
to show that the C-type orbital structure is stabilized in 
the A-AjE, magnetic phase when i| a is a small positive 
number ,Ell as it occurs in the real materials. 

In this work the Coulomb interactions (intra-orbital 
Coulomb U, inter-orbital Coulomb U', and inter-orbital 
exchange J) have been neglected. For Mn-oxides, they. 
are estimated as U — 7eV, J = 2eV, and U' = 5eV,H 
which are large compared to t. However, the result for 
the optimized distortion described in this paper, obtained 
without the Coulomb interactions, is not expected to 
change since the energy gain due to the JT-distortion 
is maximized when a single e g -electron is present per 
site. This is essentially the same effect as produced by 
a short-range repulsion. In fact, it has been checked ex- 
plicitly by using the Exact Diagonalization method that 
the JT- and breathing-distortions are not changed by V 
on a 2 x 2 cluster using the F-state in which U and J 
can be neglected. In addition, the MC simulations show 
that the probability of doble occupancy of a single or- 
bital is negligible where the A-type spin, C-type orbital 
state is stable£3 Based on all these observations, it is 
believed that that the effect of the Coulomb interaction 
is not crucial for the appearance of the A-AF state with 
the proper orbital order. Another way to rationalize this 
result is that the integration of the JT-phonons at large 
A will likely induce Coulombic interactions dynamically. 

Finally, let us briefly discuss transitions induced by the 
application of external magnetic fields on undoped man- 
ganites. When the A-AF state is stabilized, the energy 
difference (per site) obtained in our study between the 
A-AF and F states is about t/100. As a consequence, 
magnetic fields of 20 ~ 50T could drive the transition 
from A-AF to F order accompanied by a change of or- 
bital structure, interesting effect which may be observed 
in present magnetic field facilities. 

Summarizing, using numerical techniques at n = 1 it 
has been shown that the A-AF state is stable in mod- 
els with JT-phonons, using coupling values physically 
reasonable for LaMnOg. Our results indicate that it is 
not necessary to include large Coulombic interactions in 
the calculations to reproduce experimental results for the 
manganites. Considering the small but important effect 
of the octahedra tilting of the real materials, the C-type 
orbital structure (with the alternation pattern of 3x 2 — r 2 
and 3y 2 — r 2 orbitals) has been successfully reproduced 
for the A-AF phase in this context. 
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